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As a method beyond the mean-field analysis, a matrix product state (MPS) with incommen- 
surate periodicity is applied to detect phase transitions accompanied with periodicity change, 
where the incommensurate MPS is generated by acting local-spin-rotation operators with the 
incommensurate periodicity on a uniform MPS. As a commensurate/commensurate change, 
we calculate the partial ferro ~ perfect ferro phase transition in the S* = 1/2 Heisenberg model 
and its critical exponent of the magnetization curve. As a commensurate/incommensurate 
change, we calculate the S = 1 Heisenberg model with bilinear and biquadratic interactions 
which has periodicity change in the spin-spin correlation function. 
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Introduction To detect the quantum phase transition and to characterize the ground 
state of the quantum many-body Hamiltonian are highly challenging tasks. The mean-field 
analysis based on the one-body approximation is one of successful and useful methods but 
completely neglects quantum entanglement between two particles, which is quite important 
to represent the spin singlet state. The matrix product state (MPS) [1, 2] is one way to 
handle quantum entanglement gradually starting from the mean-field analysis. Especially, in 
one-dimensional (ID) quantum spin systems, many numerical studies by use of the density 
matrix renormalization group (DMRG) method [3-5] reveal that the MPS becomes a good 
variational state of finitely correlated systems. An MPS in spin 5 = 1/2 L-site systems 
is represented as = Z^o- Tr [J^^-^-^ Aj-'] |cr), where aj =t;i) f = {o"i---o'l} and A"-^ 
is an m X m matrix at j-th site. This MPS is called the vertex-type MPS. [6, 7] One of 
the variants of the MPS is called the interaction-round-a-face-type (IRF-type) MPS [6,7]: 
l^i) = Tr [nj=i ^J^''^-'^^] |<t), where cfl+i = ci and A^J'"^^^ is also an m x m matrix at 
j-th site. The vertex-type MPS is a special case of the IRF-type MPS, namely = ^J-' . 

Hereafter, we will focus on the IRF-type MPS. The MPS can express a translation-invariant 
state in the thermodynamic limit (L — > oo) with finite degrees of freedom using a uniform 
matrix j^.^''^^^^ = j^^j^^i+T-^ in spite of the infinite dimension of the Hilbert space in the 
thermodynamic limit. The infinite time-evolving block decimation (iTEBD) [8] uses this 
advantage, where the Suzuki- Trotter decomposition is used for optimization of the MPS and 
breaks the translation-invariance of the MPS. [8,9]. In order to discuss the translational 
period precisely, the Suzuki- Trotter decomposition is not used in this work. 

Even if the translation-invariance of the exact ground state in a spatially-homogeneous 
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Hamiltonian is broken spontaneously, the uniform MPS may keep translation-invariance in 
the hmit m ^ oo using the hnear combination of the degenerated ground states, as checked 
for an anti-ferromagnetic Ising model [10]. However, for fixed and finite m, the optimum 
variational state to represent g-site periodic state is the g-site modulated MPS with J^J''^^^^ = 
^modjj^q] I-*^^]' "wlicrc the commensurate period g is a positive integer and the MPS with the 
g = 1 period means the uniform, i.e., translation-invariant, MPS. Since the degrees of freedom 
of the g-site modulated MPS is g-times larger than the uniform MPS, ground states with long 
period g require large computational memory to be represented by MPSs. In this meaning, 
incommensurate period is impossible to be represented by the finite degrees of freedom. 

To overcome this problem, we have introduced the MPS with incommensurate period 
by acting local-spin-rotation operators recently [11]: |*i((5, n)) = [ W- exp(— ■jQjSj • n)] 
where a parameter Q and a unit vector n mean a pitch of the rotation angle and a rotational 
axis in the spin space, respectively. The operator Sj represents the spin operator at j-th 
site. Note that a uniform MPS is a special case of the incommensurate MPS (IC-MPS) 
|^i(Q, n)), namely = |^i(0, n)). We have already confirmed that this IC-MPS can detect 
the incommensurate properties in spin-spin correlation functions in S" = 1/2 zigzag chain. [11] 
However, it is not well understood how the transition point and the critical phenomena, for 
example critical exponents, are captured by the IC-MPS. 

In this work, we apply the IRF-type IC-MPS to following two problems about the critical 
phenomena at the phase boundary. One is (i) the commensurate/commensurate change at the 
phase transition between partial and perfect ferro-magnetism in an 5 = 1/2 antiferromagnetic 
(AF) Heisenberg chain. The other is (ii) the commensurate/IC (C/IC) change in an 5 = 1 
spin system with bilinear and biquadratic interactions. The mean-field analysis must be 
reproduced by our method in the case of m = 1 and vertex-type IC-MPS. The rigorous result 
will be obtained in the limit m — >■ cxd. Our motivation is to interpolate the two limits. 
Optimization Condition Let us show details of our numerical calculations before the 
discussion of our results. In this work, we use the IRF-typc MPS |'I'i(Q,n)) with = 
^<^j>^j+i. The m X m complex matrix-elements, pitch Q, and axis n are optimized so that 
the variational energy (^'i((5, n)|^|*i((5, n))/(*i(Q, n)|*i((5, n)) for the Hamiltonian H be- 
comes minimum by using the modified Powell method [12]. Additionally, 10-2000 initial dif- 
ferent matrices are prepared and optimized in each Hamiltonian parameter to avoid obtaining 
a false variational energy. The rotational axis n is taken as (0, 0, 1) in following results. 
Commensurate/Commensurate Change First, we evaluate a critical exponent of the 
magnetization near the phase transition between partial ferro - perfect ferro in the S = 1/2 
AF Heisenberg chain under the uniform longitudinal magnetic field: Hi = X]^ (sj-Sj_|_i — 
where Sj and mean the spin 5 = 1/2 operator at j-th site and the longitudinal magnetic 
field, respectively. The rigorous result [13] shows that the transition point is at = 2 with 
the exponent f3 = 1/2, while the mean-field calculation shows (3=1. 

As a numerical result, the magnetization curve of the 5 = 1/2 Heisenberg chain by using 
the IC-MPS is shown in Fig. 1. The magnetization per site is represented by M. The pitch 
Q becomes vr after the optimization irrespective of the magnetic field . [11] The broken 
line is calculated by the mean-filed analysis, which is M = This mean-filed solution 

is also obtained if we use the vertex-type IC-MPS with m = 1. In contrast, the IRF-type 
IC-MPS with m = 1 in this figure is better than the mean-field clearly, because the IRF-type 
MPS is able to generate non-zero entanglement entropy even if m = 1. However, the IC-MPS 
with m = 1 overestimates the magnetization in all partial-ferro region. The overestimation 
is remarkably suppressed with increasing m as shown in Fig. 1. In the case of m = 3, the 
estimation error in the magnetization becomes less than one percent. 

We evaluate the critical exponent (5 using the following relation: 1— M/Mg = — 
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Fig. 1. (Color online) Magnetization curve 
as a function of for Hi. 
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Fig. 2. (Color online) Asymptotic behavior 
of magnetization curve approaching the satu- 
rated magnetization Ms with the double log- 
arithmic plot for Hi. 



Figure 2 shows the asymptotic behavior of the magnetization approaching the saturated 
magnetization Ms = 1/2. The exponent /3 evaluated by the IC-MPS with m = 1 agrees with 
that from the mean-field analysis, namely /3 = 1, within error of five percent. The exponents 
f3 evaluated by the IC-MPS with m = 2 and 3 consist with that from the rigorous analysis, 
namely /3 = 1/2, within error of five percent for the region depicted in Fig. 2. It can be 
concluded that the IC-MPS with m = 2 is sufficiently accurate for the present model while 
the IC-MPS with m = 1 shows similar behavior to the mean-field result. 

The previous study on the Ising chain with a uniform transverse magnetic field [9] suggests 
that the vertex-type MPS leads to mean-filed behavior for any finite m in the vicinity of the 
transition point because the MPS can represent only finitely correlated state with a finite 
m. However, we cannot find such a crossover behavior for the IC-MPS with m = 2 of the 
present model in the asymptotic behavior of 1 — M/M^ down to 1 — /i^//i^ ~ 1 x 10~^. 
Commensurate/Incommensurate Change To apply our method to the C/IC change, 
we consider the S = 1 spin chain with bilinear and biquadratic interactions: H2 = cosOSj- 

Sj+i -|-sin 9{Sj ■ Sj+i)^, where Sj is the spin 5 = 1 operator at j-th site, and 9 determines the 
ratio of intensity between bilinear and biquadratic interactions. Previous numerical works 
reported that the IC pitch of spin-spin correlation functions approaches a commensurate 
pitch with increasing 6* up to 0c = 7r/4. [14, 15] This point of 9c is also a gapfull-gapless 
phase transition point, and k = ib27r/3 spin quadrupolar correlations are dominant m. 6 = 
[7r/4,7r/2]. 

As a numerical result, the change of periodicity Q as a function of the model parameter 6 
is shown in Fig. 3. The broken line represents the pitch obtained from the classical vector spin 
model, which is 2cos(5 = — cotan^. The commensurate-correlation region cannot be detected 
by the classical vector analysis. On the other hand, in the analysis with the IC-MPS, the 
C/IC change is revealed immediately with increasing m. The pitch Q/vr in the commensurate 
region agrees with 2/3 within error of 0.5 % for m = 3. Thus, the pitch Q detects the property 
change of the ground-state's correlations in this model. It can be concluded that our method 
detects the C/IC change with increasing m. 

Summary We apply a variational method with the IRF-type IC-MPS to phase transitions 
in (i) the 5 = 1/2 AF Heisenberg chain under magnetic field and (ii) the 5 = 1 AF Heisenberg 
chain with bilinear and biquadratic interactions. As the result, our method has successfully 
identified the exact phase transition points for both models. Since the MPS is identical to 
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Fig. 3. (Color online) Incommensurate pitch Q as a function of the model parameter 9 for H2- 
Characters C and IC are abbreviations of commensurate and incommensurate, respectively. 



the mean-field analysis in the simple limit and will become the exact ground state in the 
limit m — >■ oo, this method has interpolated the tvi^o limits successfully as shown in Fig. 1. 
Especially, we find that the exact critical exponent (3 = 1/2 is obtained with small m. This 
becomes a useful property of this method. Moreover, in the bilinear-biquadratic Hamiltonian, 
we have succeeded to detect the C/IC change by small m, while the analysis of the classical 
vector spin model fails to detect it. The pitch Q of the commensurate region is immediately 
converged with increase of m to Q = 27r/3, which is consistent with the period of the dominant 
spin-spin correlations in the previous studies. In the incommensurate region, the convergence 
of Q with respect to m is not enough. The improvement of optimization method to treat 
large m is one of the remaining issues. In addition, application to the transverse filed Ising 
model is one of the future issues in order to compare the previous study. [9] 
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